Na2lrC>3 as a molecular orbital crystal 
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Contrary to previous studies that classify Na2ir03 as a realization of the Heisenberg-Kitaev model 
with dominant spin-orbit coupling, we show that this system represents a highly unusual case in 
which the electronic structure is dominated by the formation of quasi- molecular orbitals (QMOs), 
with substantial quenching of the orbital moments. The QMOs consist of six atomic orbitals on an 
Ir hexagon, but each Ir atom belongs to three different QMOs. The concept of such QMOs in solids 
invokes very different physics compared to the models considered previously. Employing density 
functional theory calculations and model considerations we find that both the insulating behavior 
and the experimentally observed zigzag antiferromagnetism in Na2ir03 naturally follow from the 
QMO model. 

PACS numbers: 75.10.-b,75.10. Jm,71.70.Ej,71.15.Mb 
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High interest in the recently synthesized hexagonal iri- 
dates is due to the hypothesis 0, [H[ that the elec- 
tronic structure in these materials is dominated by the 
spin-orbit (SO) interaction. In this case, the Ir ti g bands 
are most naturally described by relativistic atomic or- 
bitals with the effective angular moment, j e fj = 3/2 and 
jeff = 1/2. In this approximation, the splitting between 
the 3/2 and 1/2 states is larger than their dispersion. The 
upper band j e ff = 1/2 is half filled and Ir atoms can be 
described as localized (j e ff = 1/2, M = 1 /is) magnetic 
moments [6] with the exchange interaction strongly af- 
fected by SO coupling. In particular, this picture leads to 
a very appealing framework known as Heisenberg-Kitaev 
model @, [![, with highly nontrivial physical properties. 
However, experimental evidence for the j e ff scenario is 
lacking 0. 

In this Letter, based on ab initio density functional 
theory (DFT) calculations and model considerations, we 
show that this picture does not apply to the actual 
Na2lr03. Instead, this system represents a highly un- 
usual case where the formation of electronic structure is 
dominated by quasi- molecular orbitals (QMOs), which 
involve six Ir atoms arranged in a hexagon. What dis- 
tinguishes this picture from molecular solids is that there 
is no associated spatial clusterization, but each Ir atom 
(via its three t^g orbitals) participates in three different 
QMOs, yet in the first approximation there is no inter- 
QMO hopping the thus formed bands are dispersionless. 

Such an electronic structure calls for a new approach. 
There is no known recipe for handling its magnetic 
properties, or adding Coulomb correlations, for instance. 
While we will not present a complete theory of spin dy- 
namics and correlations in the QMO framework, we will 
outline the general directions and most important ques- 
tions, in the expectation that this will stimulate more 
theoretical and experimental work and eventually gen- 
erate more insight. Yet, the key observable features of 



Na2lr03: small magnetic moment, unusual zigzag anti- 
ferromagnetism, and Mott-enhanced insulating behavior, 
are naturally consistent with the QMO framework. 

The main crystallographic element of Na2lr03 (see SI) 
is an Ir 4+ (5d 5 ) honeycomb layer with a Na 1+ ion located 
at its center. Each Ir is surrounded by an O octahedron, 
squeezed along the cubic [111] (hexagonal z) axis. There- 
fore, Ir d- states are split into an upper e g doublet and a 
lower t2 g triplet. The [111] squeezing further splits the ti g 
levels into a doublet and singlet; initially this effect was 
neglected 0,0,11], however, it was later included lol. [ll] 
(and overestimated) to explain the observed deviations 
from the Heisenberg-Kitaev model. 

In the previous works, after identifying the ti g — e g 
splitting it was assumed that the energy scales are W 
< (Jh, A) < E/, where W ~ 4t is the d-electron band 
width, t the effective hopping parameter, Jh the Hund's 
rule coupling, A the SO parameter, and U the on-site 
Coulomb repulsion. In this limit, the electrons are lo- 
calized and the system is a Mott insulator. While A ~ 
0.4-0.5 eV for 5d ions, the bandwidth for 5d orbitals is 
1.5-2 eV and U ~ 1 — 2 eV, Jh ~ 0.5 eV, reduced com- 
pared to typical U ~ 3 — 5 eV and Jh ~ 0.8 — 0.9 eV for 
3d electrons. Many-body renormalization may narrow 
the bands by a factor (m*/m); however, given that in Ir 
U ~ W, it is unrealistic to expect a large renormalization. 
Therefore, the usual starting point W < (J^,A) < U is 
not valid here, rather, the system is close to an itinerant 
regime, i.e., one cannot justify reducing the description 
of Na2lr03 (and possibly other iridates) to an effective 
j = 1/2 model, decoupled from the other j e ff states. 

Thus, the first step (usually skipped) is to understand 
the non-relativistic band structure. We have therefore 
performed DFT calculations (see SI) initially without SO 
effects (see Fig.[lJ solid purple lines). Inverting the band 
structure results (see SI), we obtained the corresponding 
tight-binding Hamiltonian. The leading channel (by far) 
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FIG. 1: (color online) Electronic structure of the non- 
magnetic Na2ir03 for the experimentally determined [l2[ 
crystal structure. The calculations were performed with the 
full potential local orbital (FPLO) basis using the generalized 
gradient approximation (see SI). The solid purple and dotted 
green lines refer to calculations without and with SO interac- 
tion, respectively. Note that the Fermi levels (shown by the 
horizontal dotted lines) are not aligned. 



is the nearest neighbor (NN) O-assisted hopping between 
unlike orbitals (see Fig. [2]). This was also correctly iden- 
tified previously [3, [J . There are three different types 
of NN Ir-Ir bonds; for one (we name it xy bond) (see 
Fig. 3) this hopping is only allowed between d xz and d yz 
orbitals, for the next (xz) between d yz and d xy orbitals 
and for the third bond (yz) between d xy and d xz . In our 
calculations this hopping, t[ (the prime indicates that 
the hopping is via O) is about 270 meV. Perturbatively, 



this term is proportional to t 2 d7 J(Et 2g — E p ), where p 
stands for the O p states. Ref. [5[ pointed out another 
(next nearest neighbors, NNN) O-assisted term, which 
we find to be ~ 75 meV. Jackeli and Khalliulin [4| in- 
voked another NN hopping process, between like orbitals 
pointing directly to each other. Despite the short Ir-Ir 
distance, these matrix elements are surprisingly small, 
< 30 meV. Finally, some authors lQlE] addressed the 
trigonal squeeze, which creates non-zero matrix elements 
between the same-site ti g orbitals. 

The main feature of the calculated non-relativistic 
band structure (see Fig. [1]) is formation of a singly de- 
generate (not counting spins) band state at ~ —1.2 eV, 
a doubly degenerate one at —0.7 eV, and a three-band 
manifold between —0.3 and 0.2 eV. This clear separation, 
of the order of 0.3 eV, cannot be related to the trigonal 
squeeze, as this can only split the 6 t^g bands (there are 
two Ir per cell) into a doublet and quartet. 

In order to understand this, we start with the domi- 
nant hopping, the NN O-assisted t[. Let us consider an 
electron on a given Ir site in a particular orbital state, 
say, d xz . The site has three NN neighbors. As discussed 
above, this electron can hop, with the amplitude to a 



FIG. 2: (color online) Most relevant O p- assisted hopping 
paths in idealized Na2lr03 structure. For each of the three Ir- 
Ir bond types only hopping between two particular ti g orbitals 
is possible. The same holds for the second and third nearest 
neighbor hopping via O p and Na s orbitals. Ir-Ir bonds are 
color coded as follows: xy bonds are shown by blue lines, xz 
bonds by green, and yz bonds by red ones. 




neighboring state of d yz symmetry, located at a partic- 
ular NN site. From there, it can hop further into a d xy 
state on the next site, and so on (see Figs. [2] and [3]). At 
each site, the electron has only one bond along which it 
can hop. Following the electron around, we see that after 
six hops it returns to the same state and site from where 
it started. This means that in the NN t[ approximation 
every electron is fully localized within 6 sites forming a 
hexagon. Such a state could be called a molecular orbital, 
except that there are no spatially separated molecules on 
which electrons are localized. Each Ir belongs to three 
hexagons, and each Ir-Ir bond to two. Thus, three differ- 
ent t2 g orbitals on each Ir site belong to three different 
"quasi-molecular" orbitals (QMO) and these QMOs are 
fully localized in this approximation (Fig. [3]). 

Six QMOs localized on a particular hexagon form six 
levels, listed in Table HI grouped into the lowest B in sin- 
glet, the highest Ai g singlet, and two doublets Ei g and 
E2 n . The energy separation between the lowest and the 
highest level is 4^, which is close to the calculated total 
non-relativistic t^g band width. 

We now add the O-assisted NNN hopping t f 2 . Here 
there are several such paths. However, the dominant 
hopping takes advantage of the diffuse Na s orbital (see 
Fig. ED, and is proportional to t 2 pd7r t 2 sp /(E t2g -E p ) 2 (E t2g - 
E s ) < 0. It connects unlike NNN t2 g orbitals that belong 
to the same QMO, and therefore retains the complete lo- 
calization of individual QMOs. It does shift the energy 
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FIG. 3: (color online) (a) Schematic plot of a IreNa hexagon. 
We use the same color coding as in Fig. [2j xy bonds are shown 
by blue lines and d xy orbitals by blue dots, etc. (b) A quasi- 
molecular composite orbital on a given hexagon, (c) Three 
neighboring quasi-molecular orbitals. 



TABLE I: Six quasi- molecular orbitals formed by the six tig 
atomic orbitals on a hexagon, (uj — exp(i7r/3)). Note that 
t[> and t' 2 <0 
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levels though, as shown in Table HI The upper singlet and 
doublet get closer and the lower bands move apart pro- 
viding the average energy separations of ~ 0.5, ~ 0.6, and 
~ 0.1 eV among the calculated non-relativistic subbands 
(at l^i/^l = ^ the upper two levels merge; in reality, 
l^i/^l ~ 3.3). Given that the subband widths are 0.2- 
0.3 eV, obviously, the upper doublet and singlet merge 
to form one three-band manifold. 

Several effects contribute to the residual dispersion 
of the QMO subbands. The trigonal splitting plays a 
role, albeit smaller than often assumed: the trigonal hy- 
bridization is A « 25 meV (the splitting being 3A). This 
may seem surprising, given the large distortion of the O 
octahedra, however, in triangular layers several factors of 
different signs contribute to A, and strong cancellations 
are not uncommon [l3|. Trigonal splitting, combined with 
various NN and NNN hoppings not accounted for above, 
all of them on the order of 20 meV, trigger subband dis- 
persions of 200-300 meV (see SI for further discussion). 

We shall now address the SO interaction. The corre- 
sponding bands and density of states (DOS) are shown 



in Fig. [U The lowest two subbands hardly exhibit any 
SO effect, even though the spin-orbit parameter A in Ir 
is ~ 0.4-0.5 eV, larger than both the subband widths 
and subband separation. However, a simple calculation 
shows that not only are the orbital momentum matrix 
elements between the QMOs on the same hexagon zero 
(this follows from the quenching of the orbital momentum 
in the QMO states), but they also vanish between the 
like QMOs, located at the neighboring hexagons, such 
as Biu — B\ u . Furthermore, at Y the matrix elements 
between the two lowest subbands, B\ u and E\ g , van- 
ish because of different parities; away from the Y point 
the effect of SO increases, in the first approximation, as 
F(k) = sin 2 kA + sin 2 kB + sin 2 kC, where A, B and 
C are the three vectors connecting the centers of the 
hexagons, as can be worked out by applying the L • S 
operator to the corresponding QMOs. 

The situation becomes more complex in the upper 
manifold, where three bands, A\ g and two Em, come very 
close. Even though the diagonal matrix elements, as well 
as nondiagonal elements at Y still vanish, the fact that 
A\ g and E2 U are nearly degenerate in energy induces a 
considerable SO effect at all other k-points (which grows 
linearly with /c, as y/F(k)). Note that deviations from the 
minimal model (ti,^) an d SO coupling with the lower 
Eig states also affect the bands at k = 0. We also re- 
mind that the orbital moment of the individual electronic 
states can only be finite if the QMOs mix (which is the 
case), and the direction of the orbital moment is differ- 
ent in different parts of the Brillouin zone: along one of 
the three cardinal in-plane directions it is parallel to the 
cubic x, along another to y, etc. Since the spin moment 
tends to be parallel to the orbital moment, SO is com- 
peting with the Hund's rule coupling and suppresses the 
tendency to magnetism. 

Let us now discuss the effect of the Hubbard correla- 
tions. It was initially conjectured that Na2lr03 was a 
Mott insulator. This seems counterintuitive, since sim- 
ilar Ad Ru and Rh compounds are correlated metals, 
and more diffuse bd orbitals have a smaller Hubbard 
U ~ 1.5 — 2 eV and stronger hybridization. It is hard 
to justify that this U can drive a 5/6 filled band of a 
similar width into an insulating state. Recently another, 
more logical concept has gained currency: on the LDA 
level Na2lr03 is a semimetal, barely missing being a semi- 
conductor, and a small Hubbard U just helps to enhance 
the already (spin-orbit driven) existing gap. Indeed, in 
our calculations the minimal gap is —8 meV, but the av- 
erage direct gap is 150 meV, consistent with the optical 
absorption [14}. The minimal direct (optical) gap is 50 
meV, so it is plausible that it is somewhat enhanced by 
correlation effects. 

In order to include the effect of an onsite Hubbard U 
in the QMO basis, a Z7qmo ~ U/6 has to be applied to 
each QMO[15j, with a residual Coulomb repulsion be- 
tween neighboring QMOs, V QM o ~ U/1S = ^qmo/3 
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(note that two QMOs overlap on two sites). Overall, we 
expect that the effect of the Coulomb repulsion in our sys- 
tem is similar to that in a single-site two-orbital Hubbard 
model at half filling (the upper QMO band is half-filled) 
and Uqmo ~ W ~ 150 — 200 meV. In this case, since 
Uqmo does not compete with one-electron hopping any 
more, one should expect that the gap will be enhanced 
by a considerable fraction of Z7qmo, which is consistent 
with the experiment. Thus, Hubbard correlations are of 
no qualitative importance, and only moderately enhance 
the existing gap. 

Since all electrons are fully delocalized over six sites, 
any model assuming localized spins (whether Heisenberg 
or Kitaev) is difficult to justify. On the other hand, the 
QMOs are not magnetically rigid objects and neighboring 
QMOs overlap on 2 out of 6 sites, which makes a model 
with magnetic moments localized on QMOs equally un- 
suitable [161 ] . 

We will consider therefore magnetism in the itinerant 
approach. In the non-relativistic case, the non-magnetic 
DOS shows a high peak at Ep due to E2 U and A\ g merg- 
ing and rather flat band dispersion (see Fig. [I]). Such a 
system is unstable against ferromagnetism (FM) and the 
peak is easily split gaining exchange energy (1 /ie/Ir) 
with little loss of kinetic energy. The resulting FM state 
is half-metallic (Fig. E) (see SI). 




energy (eV) 

FIG. 4: Non-relativistic non-magnetic (purple) and ferromag- 
netic (orange) density of states (DOS) of Na2lr03 calculated 
with the FPLO basis. 

Turning on the SO interaction has a drastic effect on 
magnetism. SO competes with the Hund's rule that fa- 
vors all onsite orbit als to be collinear. The spin moment 
is then reduced from 1 /iB to « 0.4 /ie/Ir for ferro-, and 
« 0.2 /ifi/Ir for the zigzag and stripe antiferromagnetic 
(AFM) arrangements (see SI). The orbital moment is par- 
allel to the spin one, reminiscent of the j e ff = 1/2 state, 
and is roughly equal in magnitude and not twice larger, 
as it should be for j e ff = 1/2. The energy gain for the FM 
case drops to a few meV/Ir [I?}, and the zigzag pattern 
evolves as the most favorable AFM state. 

Qualitatively, two closely competing ground states 
emerge from the relativistic DFT calculations: ferromag- 
netic and zigzag. In the context of an itinerant picture, 



we can argue as follows. SO creates a pseudogap at the 
Fermi level in the non-magnetic calculations (see Fig. [I]). 
This gains one-electron energy and any AFM arrange- 
ment that destroys this pseudogap incurs a penalty. From 
the three considered AFM states only zigzag preserves 
(even slightly enhances) the pseudogap (see SI). That 
gives this state an immediate energetical advantage and 
leads to the energy balance described above. Two notes 
are in place: first, all the above holds in LDA+U cal- 
culations with a reasonable atomic U (we have checked 
U up to 3.8 eV). The role of U in these systems - as 
stated previously- is merely enhancing the existing SO- 
driven gap. Second, if the DOS indeed plays a decisive 
role in magnetic interactions, it is unlikely that they can 
be meaningfully mapped onto a short-range exchange 
model, Heisenberg or otherwise. 

Summarizing, our DFT calculations demonstrate that 
Na2lr03 is close to an itinerant regime. The electronic 
structure of this system is naturally described on the ba- 
sis of quasi-molecular orbitals centered each on its own 
hexagon. This makes this, and similar materials rather 
unique. Proceeding from this description one can un- 
derstand the main properties of Na2lr03, including its 
unique zigzag magnetic ordering with small magnetic mo- 
ment. 

However, the main goal of our work is not a complete 
understanding of the magnetic properties of Na2lr03. 
We realize that this understanding is still incomplete and 
that full explanation of the weak antiferromagnetism, as 
well as of the magnetic response in this compound re- 
mains a challenge. Rather, we lay out the framework in 
which this challenge has to be resolved. We demonstrate 
that both the simplified (but correct) TB model proposed 
in previous studies 0, [H[ , and full ab initio calculations 
provide a framework that is best described by the quasi- 
molecular orbitals. This is an as yet unexplored concept 
(as opposed to molecular orbitals or atomic orbitals), and 
there are many open questions about how to treat corre- 
lations, magnetic response etc. in this framework, how- 
ever, it appears to be the only way to reduce the full 
12 atomic orbitals fog or their relativitsic combinations) 
problem to a smaller subspace (3x2 = 6) QMOs. 
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SUPPLEMENTARY INFORMATION 

We performed density functional theory (DFT) cal- 
culations considering various full potential all electron 
codes, such as WIEN2k (si|, ELK js2, and FPLO jsjj| 
using the generalized gradient approximation functional 
in its PBE form [S4[, and verified that the results agree 
reasonably well among different codes. Such compari- 
son is particularly important because the codes imple- 
ment the spin-orbit coupling in slightly different ways, 
employing usually unimportant, but in principle unequal 
approximations. In the non- relativistic calculations the 
core electrons were treated fully relativistically and the 
valence electrons non-relativistically (scalar relativistic 
approximation). In the fully relativistic calculations, i.e. 
with inclusion of spin-orbit coupling, all electrons were 
treated fully relativistically. We considered the C2/m 
crystal structure as given in Ref. IS5l and shown in Fig.lSll 

FIG. SI: Crystal structure of Na2ir03 in the cubic setting. 
The hexagonal direction is along the [111] direction in this 
setting. Ir, O and Na atoms are shown as grey, magenta, 
and yellow spheres, respectively. The three inequivalent Ir-Ir 
bonds are labeled according to their cubic directions. 




We used projective Wannier functions as implemented 
in the FPLO basis [S6| to determine a tight-binding (TB) 
representation for the Ir bd bands. In Figure IS2l we show 
the DFT band structure together with the bands cor- 
responding to the Wannier representation and the TB 
bands derived from this representation. 

In Fig. [S3] we present the projective Wannier functions 
for the bd orbitals of one Ir site. The Wannier functions 
exhibit the typical shape of the bd functions at the Ir 
site. Besides, they show a clear asymmetry due to Na as 
well as tails on the O sites. 

In order to analyze the contribution to the non- 
relativistic band structure of the various tight-binding 
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FIG. S2: Non-relativistic non-magnetic band structure of 
Na2ir03 (red symbols) shown together with the Wannier 
bands (yellow) and the tight-binding bands (blue). 

hopping parameters and its relation to the quasi- 
molecular orbital (QMO) picture, we present in Fig. [S4] 
the band structure that results if we restrict the tight- 
binding Hamiltonian to first neighbors (top left), up to 
second nearest neighbors (top right), up to third nearest 
neighbors (bottom left), and without restriction (bottom 
right). One can see that already the second neighbors 
model provides a good semiquantitative description of 
the band formation. 

In the next Figure [S5] we show the tight-binding band 
structures within the QMO model. In these calculations 
we have included the on-site trigonal splitting (the top 
left panel), adding the nearest neighbors t[ hopping (top 
right), then the second nearest neighbors t' 2 hopping (bot- 
tom left) and, finally, including also the third nearest 
neighbors hopping between the like orbital, which also 
proceeds through Na and does not take an electron out of 
the corresponding QMO (bottom right). The small dis- 
persion that arises for nearest neighbors is due to devia- 
tions from the perfect octahedral environment of iridium. 
Upon inclusion of second nearest neighbors, as mentioned 
in the main text, the upper doublet and singlet merge to 
form one three-band manifold. 

In Figure [S6| we show projections of the total density 
of states of Na2lr03 onto the quasi-molecular orbitals 
specified in Table [1] of the main text. The eigenvector 
matrix 

(1 1 1 1 1 1 \ 

1 u uj 2 -1 uj 4 uj 5 
la; 5 uj 4 -1 uj 2 uj 
1 cu 2 uj 4 1 uj 2 cu 4 
I uj 4 uj 2 1 uj 4 uj 2 
\1 -11-11 -1/ 

(with uj = exp(z7r/3)) is a unitary transformation that 
rotates the atomic Ir t<i g orbitals into the QMO orbital 



space. Ei g and E2 U states are perfectly degenerate in 
the nonrelativistic case (Figure 1561 (a)). When spin-orbit 
coupling is turned on (Figure I5B1 (b)). interestingly, the 
three upper bands are no more equivalent in this sense, 
with the central band being mostly A\ g , and the other 
two mostly Em- Importantly, there is hardly any mixing 
between the lower three bands and the upper three bands, 
emphasizing the fact that the low-energy physics is nearly 
exclusively defined by the upper three QMOs, and their 
mutual interaction, whether with or without spin-orbit. 
At the same time, one can, alternatively, project the 
same bands onto the relativistic orbitals, w/ = 1/2 
and j e ff = 3/2, and, as observed before [S7], the up- 
per two bands have more j e ff = 1/2 character than 
j e ff =3/2 character, but, for instance, at the Gamma 
point, only slightly so (more at some other points). Thus, 
even though the SO effects are considerable, they are not 
strong enough to reduce the problem to a two j e f f = 1/2 
model. 

The magnetic patterns considered in our non- 
relativistic and fully relativistic calculations are shown 
in Fig. Ell 

The ferromagnetic state shows in the absence of SO 
an energy gain of nearly 80 meV per Ir with respect 
to the non-magnetic solution and about half this value 
against competing antiferromagnetic states (zigzag and 
stripy phases); the simple Neel state is much higher in 
energy. Inclusion of SO changes the energetics consider- 
ably, as described in the main text, with the zigzag an- 
tiferromagnetic ordering becoming competitive with the 
ferromagnetic one, and lower in energy than the stripy 
phase. We deliberately do not discuss the calculated en- 
ergies in detail, because the energy differences involved 
are on the order of one meV per atom, which is beyond 
the accuracy of the density functional theory itself, and 
on the border of the technical accuracy of existing band 
structure codes. 

In Fig.[S8]we show the density of states for some mag- 
netic orderings considered in our fully relativistic calcu- 
lations. Note that the zigzag ordering preserves the non- 
magnetic pseudogap at the Fermi level, while the stripy 
ordering destroys it. 

Finally some considerations about the Hubbard U are 
at place. In fact, there are two ways of defining U in this 
case. As usually, the actual value of U depends on which 
orbitals it is being applied to. For instance, it is well 
known that in Fe pnictides the appropriate value of U 
acting on the Wannier functions combining Fe d and As 
p states is more than twice smaller that that acting on ac- 
tual atomic d orbitals since the screening effects change 
depending on the basis of active states considered. In 
molecular solids, such as fuller ides, the atomic value of 
U often appears completely irrelevant, and the physically 
meaningful value of U is the (much smaller) energy of 
Coulomb repulsion of two electrons placed on two molec- 
ular orbitals. In the case of Na2lr03 one has a choice of 
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FIG. S3: Projective Wannier functions for five of the ten Ir 5d bands, together with a structure showing the perspective. 



using an atomic U ~ 1.5-2 eV, realizing that the results 
will be strongly affected by the fact that electrons are lo- 
calized not on individual ions, but on individual QMOs, 
or of constructing U in the QMO basis. The former way 
is readily available in such formalisms as LDA+U but 
it may be a poor choice for the description of a system 
based on quasi-molecular orbitals. 
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FIG. S4: Band structure of Na2ir03 (red symbols) shown together with the tight-binding models that include only nearest 
neighbors (top left), up to next nearest neighbors (top right), up to third nearest neighbors (bottom left) and neighbors up to 
16 A (bottom right). 
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FIG. S5: Band structure of Na2ir03 (red symbols) shown together with the tight-binding models that involve only parameters 
compatible with the quasi-molecular orbitals. Only on-site parameters (top left), up to nearest neighbors (top right), up to 
second nearest neighbors (bottom left) and up to third nearest neighbors (bottom right). 
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FIG. S6: Density of states of Na2lr03 projected onto the six 
quasi- molecular orbitals given in Table [1] of the main text 
for (a) a nonrelativistic and (b) a relativistic calculation. 




FIG. S7: Three antiferromagnetic patterns considered in this 
paper: (a) zigzag, (b) stripy, and (c) Neel. 
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FIG. S8: Density of states, spin-orbit included, for two com- 
peting magnetic patterns compared with that for the non- 
magnetic state. 



